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Actuator saturation is one of the major issues of flight control in the high angle- 
of-attack region. This paper presents a saturation control scheme for linear parameter- 
varying (LPV) systems from an antiwindup control perspective. The proposed control 
approach is advantageous from the implementation standpoint because it can be thought 
of as an augmented control algorithm to the existing control system. Moreover, the syn- 
thesis condition for an antiwindup compensator is formulated as a linear matrix inequality 
(LMI) optimization problem and can be solved efficiently. We have applied the LPV anti- 
windup controller to an F-16 longitudinal autopilot control system design and compared 
it with the thrust vectoring control scheme. The nonlinear simulations show that an LPV 
antiwindup controller improves flight quality and offers advantages over thrust vectoring 
in a high angle-of-attack region. 


I. Introduction 

T HE flight control system of a tactical aircraft 
has different performance goals for low angle-of- 
attack and high angle-of-attack regions. For example, 
pilots desire fast and accurate responses for maneuver 
and attitude tracking in a low angle of attack scenario. 
While in a high angle-of-attack region, the flight con- 
trol focuses on the maintainability of aircraft stability 
with acceptable flying qualities. The potential of high 
angle-of-attack flight presents many challenges to the 
control designers. Because of aerodynamic surface sat- 
uration and control surface limitation, unconventional 
actuators such as thrust vectoring are suggested for 
aircraft maneuvering at and beyond the stall angle-of- 
attack. However, incorporation of additional thrust 
vectoring hardware could complicate the design of 
flight control laws in the poststall regime. 1 Robust 
multivariable control methods have been recently ap- 
plied to a variety of aircraft models 2 to demonstrate 
their abilities to fly at high angles-of-attack with the 
help of thrust vectoring control (see Ref. 3 and refer- 
ences therein). Besides control law design, another 
major issue of high angle-of-attack flight is control 
saturation. It is well recognized that actuator satu- 
ration degrades the performance of the flight control 
system and can even lead to instability. The destabi- 
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lizing effects of actuator saturation have been cited 
as contributing factors in several mishaps involving 
high-performance aircraft. 4 For this reason, various 
methods of preventing instability because of saturation 
have been examined, which include allocating control 
effectors and command scaling and prioritization. 

Antiwindup control method is a popular approach 
to control saturation that employs a two-step design 
procedure. The main idea of antiwindup control is to 
design the linear controller by ignoring the saturation 
nonlinearities first and then add antiwindup compen- 
sation to minimize adverse effects of the saturation 
on closed-loop performance. Desirable design require- 
ments for antiwindup compensation subject to actua- 
tor saturation are the closed-loop system stability, re- 
covery of the linear design specifications in the absence 
of saturation (linear performance recovery), and the 
smooth degradation of the linear performance in the 
presence of saturation (graceful performance degrada- 
tion). Like other saturation control techniques, the 
antiwindup compensator design often assumes a linear 
time-invariant (LTI) plant and models the saturation 
block as a sector-bounded nonlinearity. Then absolute 
stability conditions (such as Popov, circle theorems) 
are applied for the stability and performance analy- 
sis. 5 Their extension to nonlinear systems has already 
been developed. 

A general framework that unifies a large class of 
existing antiwindup control schemes in terms of two 
matrix parameters was proposed in Ref. 6. This 
framework is useful for understanding different an- 
tiwindup control schemes and motivates the devel- 
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opment of systematic procedures for designing anti- 
windup controllers that provide guaranteed stability 
and performance. Early results in antiwindup control 
often have the drawback of lacking rigorous stability 
analysis and clear exposition of performance objec- 
tives. Using an extended circle criterion, the synthesis 
condition of static antiwindup controllers is formu- 
lated as a linear matrix inequality (LMI) problem. 7 
A recent study in Ref. 8 has further revealed that an- 
tiwindup control for stable open-loop LTI systems can 
be solved globally as an LMI problem with the order 
of antiwindup compensator no more than the plant’s 
order. Alternatively, the Popov stability condition has 
also been applied to the antiwindup compensator de- 
sign problem. 9 However, the synthesis condition of 
the saturation controller is given in coupled Riccati 
equations, which are difficult to solve for the opti- 
mal solution. Most previous anti windup compensator 
designs are only applicable to open-loop stable LTI 
systems, limiting their usefulness for practical prob- 
lems. When the system is nonlinear and open-loop 
unstable, the control synthesis problem becomes very 
difficult to solve; therefore, global stabilization cannot 
be achieved . 10,11 However, in many control systems 
including flight control systems, the system dynam- 
ics are inherently nonlinear and their linearizations at 
some operation points are strictly unstable. 

The motivation for this research is twofold. First, 
the antiwindup control scheme for LTI plants in Ref. 
12 is generalized to linear parameter-varying (LPV) 
systems. This generalization is very important be- 
cause of the relevance of LPV systems to nonlinear 
systems. In fact, the LPV model can be thought of as a 
group of local descriptions of nonlinear dynamics. The 
antiwindup compensation augments existing control 
systems by maintaining stability and recovering con- 
trol performance when actuators become saturated. 
Second, saturation control for aircraft under large ma- 
neuver operations is critical because of safety concerns. 
The proposed antiwindup compensation can augment 
existing flight control algorithms, yielding enhanced 
reliability and an expanded flight envelope. In partic- 
ular, using an F-16 longitudinal dynamic model, we 
demonstrate good flight control performance of LPV 
antiwindup control in a high angle-of-attack scenario. 
Note that the study in this paper is focused on ac- 
tuator magnitude saturations, and it is the first, but 
an important step towards the actual application of 
antiwindup to high-performance flight control. The 
results can be easily generalized to actuator rate sat- 
urations, 13 which are more realistic for an advanced 
tactical aircraft. 

The notation in this paper is standard. R stands for 
the set of real numbers and R + for the non-negative 
real numbers. R mx ” is the set of real m x n ma- 
trices. The transpose of a real matrix M is denoted 
by M t . Ker (M) is used to denote the orthogonal 


complement of M. A block diagonal matrix with 
submatrices X\, X 2 , . . . , A' ; , in its diagonal is denoted 
by diag {AA, AA, . . . , X p }. We use S” xn to denote 
the real symmetric n x n matrices and S" x ” to de- 
note positive definite matrices. If M G S nxn , then 
M > 0 (M > 0) indicates that M is positive definite 
(positive semidefinite) and M < 0 (M < 0) denotes 
a negative definite (negative semidefinite) matrix. If 
a,,b G R, then sect [a, 6] denotes the conic sector de- 
fined by {( q,p ) : ( p — aq)(p — bq) < 0}. For x G R”, 
its norm is defined as ||x|| := ( x T x ) = . The space of 
square integrable functions is denoted by £o, that is, 
for any u G £ 2 , 


\u h : = 



u T (t)u(t)clt 


1 _ 

2 


is finite. 

This paper is organized as follows: Section II is de- 
voted to providing an LPV antiwindup compensator 
synthesis condition, which is a generalization of pre- 
vious results in Ref. 12. In section III, the LPV 
antiwindup compensator is applied to an F-16 longi- 
tudinal flight control system design and is compared 
with the thrust vectoring control scheme. Finally, we 
conclude the paper in section IV. 


II. LPV Antiwindup Control Synthesis 

The goal of antiwindup compensation is to modify 
nominal controllers so that if the signal from the con- 
troller is different than that which enters the plant, 
corrective feedback action is employed to reduce the 
discrepancy. Because it is impossible to provide a 
global stabilizing solution to the antiwindup control 
problem when the open-loop plant is unstable, one 
often needs to determine regional stability for satu- 
ration control and to design the controller gains in the 
guaranteed stability region. 10, 11 In Ref. 8, a sector- 
bounded input nonlinearity, sect [0,1], was considered 
for the open-loop stable plant and is not applicable to 
exponentially unstable systems. However, the derived 
performance and stability properties can be improved 
when the input nonlinearity is restricted to a smaller 
sector region. As a result, this modification leads to 
regional stability of the antiwindup compensated sys- 
tem and extends the antiwindup control technique to 
exponentially unstable open-loop systems. 

Consider an LPV plant P p described by 


Xp 


Ap{p) Bpi(p) Bp2(p) 


Xp 

e 

= 

Cpi(p) Dpu(p) Dpi 2 (p) 


d 

y _ 


C P 2(p) D p -2i{p) D P' 2‘2 ip) _ 


_cr(u)_ 


where the plant state x p G R”**. y G R” " is the 
measurement for control, and a(u) G R”“ is the satu- 
rated control input, e G R”» is the controlled output 
and d G R" d is the disturbance input. It is assumed 
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that the vector-valued parameter p evolves continu- 
ously over time and its range is limited to a compact 
subset V C R s . In addition, its time derivative is 
assumed to be bounded and satisfy the constraint 
i '_j < pi < Vi,i = 1,2, ...,s. For notational pur- 
poses, denote V = {v : £,; < v, < i'u * = 1,2,..., s}, 
where V is a given convex polytope in R s that contains 
the origin. Given the sets V and V, the parameter v- 
variation set is defined as 

T v v = {p E C 1 (R + , R s ) : p(t) EV, p(t ) EV, Vi > 0} 

So Tj, specifies the set of all allowable parameter tra- 
jectories. 

All matrix valued state-space data are continuous 
and have appropriate dimensions. For simplicity, we 
assume that 


(Al) (Ap(p),B P 2 (p),Cp 2 (p)) triple is parameter- 
dependent stabilizable and detectable for all p. 

(A2) The matrix functions [ B p2 (p ) Dp 12 {p)] and 

\Cp 2 (p) D P 2 i(p)] have full row rank for all p. 

(A3) D P 22(p ) = 0. 

The actuator nonlinearity under consideration is a 
piecewise-linear saturation 

, J in k| <c ax 

U ®' ( sgn(u,:)u™ ax \uj j > w’ nax 

for i = 1,2,..., n u . The antiwindup control structure 
is shown in Fig. 1(a). 

Following the standard antiwindup procedure, a 
nominal LPV controller K p is designed first by ignor- 
ing the input nonlinearity. Different control design 
techniques can be employed to achieve this goal. A 
systematic way to do this is through LPV control the- 
ory. 11 Because of assumption (Al), the nominal 
controller K p exists and is capable of stabilizing the 
open-loop system when no input saturation exists, and 
its design will determine the nominal performance of 
the closed-loop LPV system. We assume that such a 
controller is given by 


Xk 


Ak (p, p) 

Bk(p) 

Xk 

u 


. C k (p) 

Dkip) 

. v _ 


where Xk E R" fc . The variables v \ , Vo are the auxiliary 
inputs provided by an antiwindup compensator. They 
are used to condition the nominal controller when the 
control input is saturated. 

Our objective is to design an LPV antiwindup com- 
pensator AW P such that the adverse effect of input sat- 
urations are minimized in terms of induced £2 norm. 
The antiwindup compensator is in the form of 


Xaw 


Vi 

= 

1 ’2 



Aaw ( P- p) 
Caw (/b P) 


Baw ip) 
Caw ( p ) 



( 3 ) 



(a) LPV antiwindup control structure 


e 


11 


q 



(b) Equivalent transformation 

Fig. 1 Nonlinear saturation control diagram. 

where the state x aw E R” al "; the size of the compen- 
sator state will be determined in the sequel. Such 
antiwindup compensation schemes provide a compu- 
tationally efficient technique for “retro-fitting” exist- 
ing unconstrained controller K p to account for input 
nonlinearities, thereby eliminating controller windup 
problems for input saturated nonlinear systems. 

The LPV antiwindup control diagram in Fig. 
1(a) can be transformed to its equivalent form by 
substituting each actuator saturation with a dead- 
band nonlinearity A,; = 1 — a(ui)/ui and A = 
diag{Ai, A 2 , . . . , A„ u ) as shown in Fig. 1(b). This 
allows recasting of the compensator design problem 
into a robust LPV control paradigm. 



Fig. 2 Restricted range input saturation. 

It can be seen that the deadband uncertainty A, : re- 
sides in the conic sector [0,1]. If the maximal value 
of the control input signal u, is restricted to be less 
than (1/(1 — ki)) w™ ax with 0 < < 1, then the 

nonlinearity A,; is reduced to sect[0, ki] as shown in 
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Fig. 2. Consequently, we have A G sect[0, K\ with 
K = diag{/ci , k' 2 , ■ ■ ■ , k nu }. The resulting stability no- 
tion then becomes regional stability, which is more 
conservative than the actual stability limit. However, 
this restriction of uncertainty A is capable of extend- 
ing the antiwindup scheme to open-loop exponentially 
unstable systems. 

Let the system G p be the interconnection of the 
open-loop system P p and the nominal controller K p , 
but exclude the antiwindup compensator. Then its 
dynamic equation is 


x 

u 

e 
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A{p, p) 
Coip) 
Clip) 
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q = A u 


Bo(p) 
Doo(p) 
Dio ip) 
I 


Blip) 
A nip) 
Dii{p) 
0 


B*(p) 

Doiip) 

Diiip) 

0 


1 x 

q 

d 

Vi 

J V-2 


(4) 

(5) 


where x T = [x p x^] G R" with n = n p + n k and 
A G sect[0, K]. 


Aip, p) = 
Boip) = 
Blip) = 


l/i T B p iD fcC p i B p iCk 

BkCpi A k 


—Bpi 

0 


Biip) = 


Bpi + B p iD k Dpi i 
BkDpn 


0 Bpi 

1 0 


Co ip) — [DkCpi Ck\ 

Clip) = \Cp i + D p nDkC p i DpiiCk \ 

Dooip) = 0, D 0 iip) = DkDpn 

Doi(p) = [0 I] , Dioip) = —Dpii 

Dnip) = D p ii + D p iiDkD p n, Dnip) = [0 A p i 2 ] 


Note that the state-space data have linear dependency 
of parameter rate p under the LPV control design 
framework. 

Denote x J = \x T xj iw \ . Then the final closed-loop 
system T p = B^iGp, AW P ) is described by 


Xd 

u 



A c iip,p) B 0 , dip) B 1 , c iip) 
Co, clip, p) D 0 o,dip) D 0 i,dip) 


Xcl 

q 
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Cl, dip, p) D 10 ,d(p) D\ \.cl(P)_ 


d 

i 


q = A u (7) 


where !/)(■, •) stands for a lower linear fractional trans- 
formation . 17 The state-space data of the closed-loop 
system T p are related to the interconnected system 
G p and the antiwindup compensator AW p . Specifi- 
cally, the closed-loop state-space data depend on the 
antiwindup compensator gain in affine form. The fol- 
lowing theorem provides a synthesis condition for the 
antiwindup compensator. 


Theorem 1 (synthesis condition for LPV antiwindup 
compensator) 

Given scalars 0 < fc, < 1, i = 1,2 ,...,n u , the LPV 
open-loop system P p with a parameter- dependent sta- 
bilizing nominal controller K p , if there exists a pair 
of positive definite matrix functions Ruip) G S'^ pXn ' p , 
Sip) G S” x ” and a diagonal matrix function V ip) = 
diag {rq (p), • • • ,v nu ip)} > 0 satisfying eqns. (8)-(10) 
for all ip,v) G V x V, then there exists an n p th or- 
der LPV antiwindup compensator AW P to stabilize 
the closed-loop system exponentially and have the per- 
formance ||e||o < 7 1 1 r/ 1 1 o for all pi-) G Tfi when the 
condition \ ufi[ < (1/(1 — kf)) u™ ax , i = 1,2 ,...,n u 
holds. 

The proof can be found in the Appendix. Because 
only the ( 1 , 1 ) element of the R matrix function is con- 
strained in the LMIs ( 8 )— (10), it is always possible to 
augment matrix Rn to R in satisfying the preceding 
coupling condition. For example, one can choose 


Ruip) 


[i 

o] s-Hp) 

o' 

I 


[0 i]s-Hp)^ 

r 

[0 

i]s-Hp) 

0 

I 



The resulting R is positive definite because of condi- 
tion (10). Also, R{p) — S~ 1 ip) > 0 is satisfied for 
the selected R matrix function. The rank condition is 
trivially satisfied if one chooses n aw = n p . 

Recall that the LPV antiwindup compensator syn- 
thesis problem is originally formulated as a robust 
control problem. This usually leads to a nonconvex 
solvability condition as bilinear matrix inequalities. 
However, because of the special structure of the an- 
tiwindup controller, the resulting synthesis condition 
is convex in terms of matrix variables Ruip), Sip), 
and V(p). In fact, the solvability condition for the 
LPV antiwindup compensator is given as an infinite- 
dimensional LMI optimization problem, for which an 
efficient numerical algorithm exists to solve it approx- 
imately . 18 This can be done by parameterizing the 
matrix variables using a finite set of scalar basis func- 
tions as 


N f N g 

Rilip) = ^2fi(p)Ru,i, Sip) = ]gjip)Sj 

i= 1 j = 1 

N h 

vip) = Y hkip)V k 

k= 1 

where {/,(/<) }; V/ , , {.'/.,(/') } ; V ' P and {//,(/<) } A V '- , are 
user-specified scalar basis functions. Rnj,Sj, and V k 
are new optimization variables to be determined. After 
such a parameterization, the LPV synthesis conditions 
can be solved using a gridding method over the param- 
eter space. 
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Rnip)A p ip) + A p (p)Rii(p) - — I 

+2B p2 (p)K- 1 V(p)(I - K-^BUp) j 


J C p i{p)Rii{p) 
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| D p uip) 
-7 In d _ 


S(p)A(p, v) 4- A T {p , v)S(p) + 

Bf(p)S(p) 
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Rii(p) 
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[ln P 0] 
Sip) 


> 0 


S{p)Biip) 

“7 In d 
Dnip) 


Clip) 

D T uip) 

-jCle 


< 0 


( 9 ) 


( 10 ) 


After solving Rnip),Sip), and Vip) matrix func- 
tions, the LPV antiwindup compensator gain can be 
determined by either solving an LMI feasibility prob- 
lem or the antiwindup compensator can be constructed 
explicitly as shown in the following theorem. The ex- 
plicit construction approach is advantageous because 
it avoids possible numerical ill-conditioning when solv- 
ing the preceding feasibility LMI problem. Moreover, 
it connects the antiwindup controller directly to the 
plant and nominal controller gains. 

Theorem 2 (LPV antiwindup compensator construc- 
tion) 

Given the solutions Rnip),Sip), 7 , and Vip) = 
W~ l ip) of the LMIs (8)-(10), let M(p)N r ip) = I n - 
Rip) Sip) with M,N G S nx % and H T = [/„ 0] , 

then one n p -th order LPV antiwindup compensator 
AW P can be constructed through the following scheme: 

1. Compute a feasible D aw ip) G R"“ x such that 
eqn. (11) holds for all p G V, 

2. Compute the least square solutions of the 

linear equations (12)-(13) for B aw ip) G 
R"x™“, Cawip) £ R”“ xn * > and the matrix 

A-aw ip,p) G R" x " p is expressed as in eqn. (If), 

3. Convert the transformed antiwindup compensator 
gain to its original state-space data by eqn. (15). 

This proof can be found in the Appendix. The 
explicit antiwindup construction scheme can also be 
applied to open-loop stable systems. However, because 
the synthesis condition for an open-loop stable plant 
does not involve W (Ref. 8 ), we need to solve both the 
feasible D aw and W matrix functions at the first step. 
The remaining steps are the same by setting K = I. 

III. Saturation Control for Flight 
Dynamics 

In this section, the proposed LPV antiwindup con- 
trol synthesis technique is applied to flight dynamics. 


The system to be controlled is the longitudinal F-16 
aircraft model based on NASA Langley Research Cen- 
ter (LaRC) wind tunnel tests , 19 which is described by 
Stevens and Lewis in great detail . 20 

A. Aircraft Model 

The full nonlinear longitudinal model of an F-16 air- 
craft is given as follows: 


qSc 

V = [C xq (a) cosa + C zq ia) sin a] q 

+ — [( '.,.( 0 , S e ) cos a + C z (a,S e ) sin a] 

— g sm(# — a) + T (16) 

m 

a = {l + .^y 2 [Czqia) cos a - C xq (a) sin 0 ] j q 


+ “Tv [CL (a, S e ) cos a 
mV 


C x ict,S e ) sin a] 


q . m sma 

+ _ cos( 9_ a )_T— 


Q = 


qSc 


21 yV 


[cC mq ia) + A C zq (a)] q 


+ 


qSc 


C m (a, S e ) + — C z ia,S e ) 
c 


e = q 


(17) 


(18) 

(19) 


The aerodynamic coefficients are provided as look-up 
tables from LaRC wind tunnel tests on a scale model 
of F-16 aircraft. The data apply to the speed range 
up to about Mach 0.6, and cover a very wide range 
of angle-of-attack (—20 deg < a < 90 deg). However, 
conventional aerodynamic math models for use in air- 
craft simulation or flight control design have become 
increasingly deficient in the poststall region . 20 So, 
the investigation on the robust control of aircraft in 
the high angle-of-attack region with uncertain aero- 
dynamic coefficients is one of the challenging research 
topics. 
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(15) 


The F-16 is powered by an afterburning turbofan 
jet engine, which produces a thrust force in the x axis. 
The LaRC data include a model of the engine in which 
the thrust response is modeled with a first-order lag; 
the lag time constant is a function of the actual engine 
power level and the commanded power. The throttle 
position is related to the commanded power level. For 
convenience, the actual power level is also considered 
as a state variable in longitudinal dynamics. 

The state and input variables of the F-16 model are 
as follows: 


x = < 


V 

a 

q 

e 

pow 


true airspeed (ft/s) 
angle-of-attack (rad) 
pitch rate (rad/s) 
pitch angle (rad) 
actual power level (0-100) 


u = 



throttle position (0-1) 

elevator angle (deg) (—25d,eg-25deg) 


By slight abuse of notations, the variables just listed 


also represent perturbations from their equilibrium 
states when linearization is considered. In addition, V, 
q, and flight-path angle 7 are selected as outputs. To 
describe the nonlinear F-16 model by an LPV system, 
we first need to find the wings-level equilibrium solu- 
tions at several flight conditions in the design envelope. 
The local linear models are then obtained by lineariz- 
ing the nonlinear equations of motion at those equi- 
librium points. The flight envelope of interest covers 
aircraft speeds between 160 and 200 ft/s and angles-of- 
attack between 20 and 45 deg. These two variables are 
used as scheduling parameters in the LPV modeling of 
F-16 longitudinal dynamics. The points at which the 
nonlinear model is linearized are marked by a x sym- 
bol in Fig. 3. This group of linearized models consist 
of the LPV representation of the nonlinear F-16 lon- 
gitudinal dynamics within the chosen flight envelope. 
The corresponding dynamic pressure range is between 
20 and 50 lb/ ft 2 . This portion of flight envelope is cho- 
sen because the moderately high angle-of-attack and 
the low dynamic pressure cause aerodynamic control 
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surface saturation, which is a major concern in this re- 
search. The actuators are explicitly constrained. Most 
of the time, we found that the elevator angle saturation 
becomes a limiting factor for flight control effective- 
ness. 

210 
200 
190 

I 

% 180 
170 
160 


angle of attack a (deg) 

Fig. 3 Flight trim points. 



B. Problem Setup 

The design objective of the nominal LPV controller 
in this research is to track the flight-path angle com- 
mand with the tracking error about 1.25% of the com- 
mand in the steady-state motion. This kind of problem 
is conveniently formulated as a model-following prob- 
lem, 21, 22 where the ideal model to be followed is chosen 
to be a second-order filter based on desired flying qual- 
ities 

2 

Tideal ^ideal (20) 

ben id S~ + 2(5,1,.;,] • (dideal ' S ‘ w id e al 

The implicit model-following framework allows for di- 
rect incorporation of flying quality specifications into 
the control design. A block diagram of the system in- 
terconnection for synthesizing the nominal controller 
is shown in Fig. 4, where P p is the model set of lin- 
earized aircraft dynamics at different operating points 
p and n is 3-dimensional sensor noise. 

The weighting functions are chosen as 

W p = W n = dia.g {0.8, 0.6, 0.1} 

W„ = diag{l,10,T,_L} 

TI . 2.25 

ldeal s 2 + 2.4s + 2.25 


The throttle and elevator actuator dynamics are mod- 
eled as first-order lag filters. 


W^7 d 


5 

s + 5 ’ 


<W7 nd 


20 

s + 20 


Both the positions and the rates of control inputs are 
fed into W u to penalize the control effort. Therefore, 
the system matrix of W act is derived as follows. 


Wact = diag 


-5 

5 ' 


' -20 

20 ' 

1 

0 

5 

1 

0 

-5 

5 


-20 

20 


C. Design Results 

The nominal controller is designed by formulating 
an LPV synthesis problem, which can be solved us- 
ing either a single or parameter-dependent quadratic 
Lyapunov function over all gridding points in the two 
dimensional parameter space. 14,15 In this work, the 
single quadratic Lyapunov function (SQLF) approach 
is chosen to reduce the computation time in the nom- 
inal LPV control synthesis phase. The performance 
obtained through the SQLF approach is 7 n0 m = 8.0. 

The proposed LPV antiwindup scheme is then ap- 
plied to the control of an F-16 aircraft, and we also 
choose the SQLF approach to perform the LPV anti- 
windup synthesis. A series of LPV antiwindup com- 
pensators AW P are designed by gradually decreasing 
the ki(i = 1, 2) value from 1 to 0.99. Table 1 shows the 
corresponding performance level 7 and the trade-off 
between the sector-bound constraint and the achieved 
performance. Because the linearized model at one of 
the gridding points is unstable, the synthesis condition 
for LPV antiwindup control is infeasible for k = 1. 
Compared with the nominal performance, the worst 
performance level when k t = 1 — 10“° indicates the 
strong adverse effect of saturation nonlinearity. On 
the other hand, the antiwindup compensator almost 
recovers the nominal closed-loop performance 7 nom 
when ki = 0.99. The antiwindup compensator cor- 
responding to ki = 0.999 will be used in the following 
nonlinear simulations. The achieved performance by 
this LPV antiwindup compensator is 25.5 compared 
to the nominal performance of 8.0. Therefore the 
controlled performance will be sacrificed in order to 
achieve a guaranteed stability result. 

Table 1 kioo performance level vs. sector range 

[0, ki]. 


Sector range [0, ki] 

Hoc performance 7 

1 - 10 -5 

1251.92 

1 - 10- 4 

196.55 

0.999 

25.50 

0.99 

8.00 


D. Nonlinear Simulations 

The LPV nominal controller K p and antiwindup 
compensator AW P are tested first at one designed 
flight condition, which is a trimmed flight at V = 
160 ft/s and a = 35 deg. The magnitude limits of ac- 
tuators are enforced during the nonlinear simulation. 
A flight-path doublet input with magnitude ±2 deg 
is used to demonstrate the performance of the anti- 
windup compensator. Both the nominal controller and 
antiwindup compensator are fixed LTI controllers for 
the given flight condition derived from the LPV con- 
trollers K p and AW P . The responses of the nonlinear 
system are shown in Fig. 5, where the solid line is the 
trajectory with the antiwindup compensator, and the 
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Fig. 4 Open-loop interconnection for LPV control design. 


dashed line is the trajectory without the antiwindup 
compensator and controlled only by the nominal con- 
troller. The dotted line in Fig. 5 (a) represents the 
flight-path angle response of the ideal model. It is ob- 
served that the performance without the antiwindup 
compensator becomes worse around t = 15s because of 
elevator saturation. However, the insertion of the anti- 
windup compensator quickly overcomes the saturation 
and greatly improves the tracking performance. 

Another test condition is selected at V = 200 ft/s 
and a = 2‘2deg, which is not a designed flight condi- 
tion, so both the nominal controller and antiwindup 
compensator are not fixed LTI controllers but varying 
in time. However, the computation time will be in- 
creased greatly if the controllers are constructed online 
according to the current parameter value. To reduce 
computational cost, a relatively coarse gridding pa- 
rameter space V is used for synthesis, and then K p and 
AW P are constructed off-line at finer gridding points 
and saved as look-up tables. In our study, 10 points 
are used for synthesis as shown in Fig. 3, and 45 points 
for constructing look-up tables of K p and AW P . When 
doing simulations, the nominal controller and anti- 
windup compensator gains are calculated using linear 
interpolation at the current values of the scheduling 
parameters between the grid point solutions. 

Note that the nominal controller K p is designed for 
a set of linearized plant models at the gridding points, 
so that the output of K p is the control deviation from 
the nominal control maintaining trim condition, and 
the actual control input to the nonlinear plant should 
be the summation of both terms. When K p is parame- 
ter varying in simulations, the instant nominal control 
inputs should be determined at each time. This can 
be done by building up another look-up table of nom- 
inal control inputs at finer gridding points beforehand 
and interpolating in real time according to the current 
parameter values. 

For the flight condition just listed, a flight-path 
doublet input with magnitude ±4 deg is used. The 
simulation results has confirmed that the trajectory of 
parameters evolves about half of the parameter space. 
In this case, the system might not achieve the desired 
performance, even maintain stability by using a single 


LTI antiwindup compensator. Instead, it is crucial to 
use LPV nominal control and LPV antiwindup com- 
pensation to achieve good controlled performance. It is 
also shown that the variations of the scheduling param- 
eters do not exceed the assumed limits: \V\ < ‘20ft/ s' 2 
and |d| < 10 deg/s, which are used in the LPV control 
synthesis stage. Moreover, the sector bound constraint 
on control input u,; < 1000u“ ax is satisfied over the en- 
tire simulation time. Fig. 6 shows the aircraft response 
with an LPV antiwindup compensator, and the dotted 
line in Fig. 6(a) is the ideal flight-path angle response. 
It is observed that the LPV antiwindup compensator 
achieves the desired performance objective. However, 
when the antiwindup compensator is unused, both the 
elevator and the throttle are saturated severely and 
the system goes unstable. For clarity purpose, the re- 
sponse without LPV antiwindup compensator will not 
be shown in Fig. 6. Note that the variations of eleva- 
tor angle and throttle position in this case study might 
go beyond the typical rate saturation limits, and thus 
are not realistic for implementation. 

E. Comparison with Thrust Vectoring Control 

Unconventional control effectors like thrust vector- 
ing are another way to enhance maneuverability of 
the modern aircraft at and beyond the stall angle-of- 
attack. To compare with the saturation control scheme 
using an antiwindup compensator, a simple thrust vec- 
toring model is added to the previous nonlinear F-16 
model to provide additional longitudinal axis control 
power. For simplicity, there are only thrust compo- 
nents along x and z body axes. A detailed description 
of the thrust vectoring model can be found in Ref. 23. 

Now the input variables are elevator deflection S e 
and pitch thrust vectoring nozzle deflection S ptv , whose 
position limit is —17 deg < S ptv < 17 deg. Because 
the state variable pow is not related to either input in 
the linearized model, it is neglected and the number 
of states is reduced to four. The longitudinal linear 
model of the F-16 aircraft at a given trim point with 
throttle and elevator as controls and that with elevator 
and thrust vectoring as controls are 


xi = A\X\ + B\u\ 


-4 1 x I V B\ 



( 21 ) 
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(c) Elevator angle 


(d) Throttle position 


Fig. 5 Saturation control of F-16 dynamics with and without LTI antiwindup compensator. 


X-2 = A 2 X 2 + B 2 U 2 


A 2 x- 2 + B 2 



( 22 ) 


where the states x 2 include velocity V, angle-of- attack 
a, pitch rate q, and pitch angle 6. 

There are usually two methods to generate thrust 
vector commands, daisy-chain 21,23 and ganged control 
methods. 1,24,25 For both methods, a controller is de- 
signed first based on the generalized control, and the 
real control inputs are then generated using a control 
selector. The former approach commands thrust vec- 
toring only when the conventional aerodynamic con- 
trol surfaces are ineffective, i.e. they are unable to 
generate the necessary forces and moments required 
for commanded maneuvers. However, it is difficult to 
define the generalized control for high-order systems. 
The basic idea of the latter approach is to redefine 
the control contribution to the state dynamics equa- 
tion with different weights, and the thrust vectoring 
always works in flight control. The disadvantage of 
the ganged control scheme is that it cannot precisely 
determine the weights of different controls. 

The control allocation method used in this research 
is different from those two approaches because the gen- 
eralized control to design the controller is not used. 
The aircraft model with conventional control effec- 
tors is used in the stage of controller design and 
the designed controller is considered as the nominal 
controller K p . A control allocation function is then 
implemented in the nonlinear simulation and thrust 
vectoring is commanded only when the elevator un- 


dergoes saturation. So, the thrust vectoring control 
is not used for control design and only used for con- 
trol implementation. This control allocation method 
avoids defining the generalized control and setting the 
weights for different control effectors. 

In the linearized model with conventional control ef- 
fectors (21), the state variable paw is only related to 
the throttle, which has no relation with the four other 
states. For simplicity in describing the control allo- 
cation scheme, the first four states with only elevator 
deflection are considered in (21). So the states and A 
matrix of the two linearized models are the same. 


x = Ax + B{ u\ = Ax + B{ 5* 


x = Ax + B 2 u 2 = Ax 4- B 2 


SptV 

Sr 


where B represents the simplified input matrix, and 
5* is the elevator deflection. Ideally, we would like the 
actual control 


B 2 



to provide control effect equivalent to that of the artifi- 
cial control B*S* . Then the practical thrust vectoring 
and elevator angle can be obtained by solving the 
following optimization problem when the elevator is 
saturated. 


min Hflatta - B^u*^ 


subj. to 


'-17 + *°*? 

< 

dptv 

< 

"17 + Sp tv ~ 

. —25 + $e . 


_s e 


25 + 5® 
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time (s) 


(a) Flight-path angle 




(b) Angle-of-attack 



(c) Elevator angle 


(d) Throttle position 


Fig. 6 Nonlinear doublet response with LPV antiwindup compensator. 


where 6p tv and <5° are equilibrium inputs at trim con- 
dition. 

To compare the control effect of the antiwindup 
compensator and thrust vectoring, a ±2 deg flight path 
doublet is commanded at the flight condition V = 
160 ft/s and a = 35 deg, same as in Fig. 5. The solid 
line in Fig. 7 represents the response of the closed-loop 
system with thrust vectoring. Similarly, the dotted 
line in Fig. 7 is the ideal trajectory of the flight-path 
angle. A comparison between the solid responses in 
Figs. 5 and 7 illustrates that the antiwindup control 
achieves better tracking performance than thrust vec- 
toring control scheme. 

IV. Conclusion 

The goal of flight control in the high angle-of-attack 
region is to maintain aircraft stability with acceptable 
flight qualities. However, the saturation of the con- 
ventional aerodynamic control surfaces presents many 
challenges to the control designers, for example, actu- 
ator saturation is an important issue to near-stall and 
poststall flight conditions. Without cost/hardware 
concerns, unconventional actuators such as thrust vec- 
toring are usually suggested to compensate the conven- 
tional aerodynamic control surfaces. However, the in- 
corporation of additional thrust vectoring nozzle could 
complicate the design of flight control laws in the post- 
stall regime. The antiwindup method provides an 
alternative approach to handle control saturation. Ap- 
plying antiwindup control scheme to flight control is 
promising because no additional actuator is needed to 


compensate control authority. The implementation of 
antiwindup controllers could be done by simply mod- 
ifying flight control software. 

An antiwindup control scheme for LPV systems has 
been developed in this research. The extension of the 
antiwindup control idea to LPV systems provides a 
practical approach for nonlinear flight dynamics in the 
presence of actuator saturation. Because of the special 
structure of the antiwindup control scheme, the LPV 
antiwindup control synthesis condition is solvable by 
LMT optimizations. The LPV antiwindup compensa- 
tion has been applied to an F-16 aircraft and compared 
with the thrust vectoring control scheme. By aug- 
menting the nominal longitudinal autopilot with an 
antiwindup compensator, it has been shown through 
nonlinear simulation that the F-16 aircraft maintains 
stability and adequate control performance in case of 
actuator saturation, while the performance achieved 
by the thrust vectoring control scheme is undesirable 
when a large maneuvering operation is commanded. 


Appendix: Proof of Main Results 

Theorem 1 

Proof: Consider a Lyapunov function V(x) = 

x T X c [(p)x for the closed-loop system T, then a suf- 
ficient condition for the exponential stability and per- 
formance can be established from the inequality 


V + — e T e — - (I 1 (I + 2q T W(Ku 


q)< 0 
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time (s) 


(a) Flight-path angle 



(c) Elevator angle 



(b) Angle-of-attack 



(d) Thrust vectoring nozzle deflection 


Fig. 7 Nonlinear doublet response with thrust vectoring. 


using 5-procedure. Note that it is equivalent to 


where 


' A T d X cl +X cl A c A 

{ } 

Bl cl X cl + WK C 0 , cl 


Ci, d 

A,/ II i,, i 
WKD 0hcl 
-7 1 
Du, d 


C T ■ 

°1 ,cl 

Dl 0 , c. 
Dfi, c 

-l 1 . 


XciB 0 ,c + Cl d KW 

(WKD 00 ,d + r/> md KW\ 
\ -2 W J 


—2W 

T 

01 ,c 

Dio, d 


Dli, d KW 


< 0 


(Al) 


where the closed-loop state-space data are as follows: 


A(p,p) = 

B 0 (P) = 


Mp, p) 

0 

Bo(p) 

0 


Blip) 

Blip) 

Blip) 


[Qi Q-2 Q 3 ] — 


Co(p) = [C 0 (p) 0] , 

0 B 2 (p) 


0 


Bi (p)= B ^ P) 
Clip) = [Clip) 


0 D 0 2ip) 


0 Di^ip) 


0 

I 

0 

1 0 1 

0 

0 

I 

1 

0 


0] 


c„ 


uiP,f>) 

uiP,P) 


Baw 

D n in 


ip) 

(P)J 


U(p,p) 

Bo, dip) 

Bi,d 

ip) 1 


D, clip, p) 

D 00 , dip) 

D 01 , ci 

'ip) 


1 ,cliP,i>) 

Dio, clip) 

D 11 , a 

'ip). 


Aip,p) 

Boip) 

Blip)' 


\Blip)] 

C-oip) 

Dooip) 

Doiip) 

+ 

Blip) 

. Clip) 

Dioip) 

Dn(p)_ 


[Bl(p) \ 


[Qi Q-2 Qz\ (A2) 


Denote 0(p, p) = 


Ap,p) 

:iP,(>) 


Ba 

D n 


vip) 

vip) 


and Vip) = 


W 1 (p). Inequality (Al) can be rewritten as follows: 


yip,p) + V 1 (p)Qip,p)Q+Q 1 Q 1 ip,p)V(p)<0 (A3) 
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for all ( p , p) E V x V with 


shown that 


*(p,p) = 


A T X C 1 + -Y c /»4 + 'Ypi 

13qX c i + WKCo 
BfX d 
Ci 


X cl B 0 +C^KW 

XciBi 

cr 

WKD 00 + Dq 0 KW) 
—2 IP J 

WKDqi 

Dio 

D ( ' n KW 

- 7 1 

Dh 

Dio 

Du 

~ll 


P(p) = [Pi-Y c , PoA'W 0 P 3 ] 

Q = [Qi 0.2 0.3 0] 


A f^(p)^(p,p)M v (p) 


RnAp 4- ApRn — Yjn 


+2B p2 K~ 1 V(I - I<~ 1 )Bp 2 
C p iRn 

+2D pl2 K- 1 V(I-K- 1 )B^ 


RuCZ. 




+2B p2 K- 1 V(I-K- 1 )DpV2[ 

-7 1 

+2D pl2 K~ 1 V(I - K~ 1 )Dj 12 


Partition the matrix function X d compatibly to the 
states of an interconnected system G p and antiwindup 
compensator AW P as n = n p + nu and n aw , and let 


Xd(p) = 


KYp) = 


\ S(p ) 

N(p)~ 

N T (p ) 

? 

' R(p) 

M(p) 

M T (p ) 

? 

Rii(p) 

RMp) 

RUp ) 

R 22 (p) 

M 

t TpT 


where M(p)N T (p) = I — R(p)S(p). According to elim- 
ination lemma, 26 inequality (A3) is equivalent to 


A fQ^{p,p)jV Q 


SA + A T S + ypi?— SB 1 Cl 
OPi 

BJS - 1 - 7 J D\ \ 

Ci Du -7 


which are the same as conditions (8) and (9), respec- 
tively. 

Given the definition of matrix functions X c j and 
-Y^ 1 , the coupling condition between R and S would 
be 


A/'p(p)^(p,p)A'V(p) < 0 and A /q^(p,p)Aq < 0 

(A4) 

for all (p, p) E V x V. M-p and Mq expand the null 
spaces of matrix functions V and Q, which are 


R(p) I 
I S(p ) 


> 0 and rank[7?(p) — S 1 (p)] < 


Mv (P) = diag { X;, 1 (p) , W- 1 (p) , /, I } 

I 0 0 

0 0 0 

0 0 0 

x -k-^b; i 2 (p) -K-ip^ip) o~ 
0 0 / 

0 10 


A'o — 


' I 

0 

0 

0 " 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

1 


for all p E V. This is equivalent to the third LMI 
because only Rn is constrained. Q.E.D. 


Theorem 2 

Proof: The derivation of the antiwindup controller 
formula basically follows the procedure outlined in Ref. 
12. Define 


Through lengthy algebraic manipulations, it can be 


1 

RH ' 

„ r 5 

H 

0 

M t H\ ’ 

h 

1 

II 

CN 

NJ 

0 
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Then it can be shown that X c iZ\ = Z 2 . Also we have 
the following congruent transformation: 


ZiXaZi = 
Zi 


T dX c i 
1 ~dF 


H t {R 


dS 


S 

h t 


dS 

dt 


H 

H t RH 


(R.j | + M— 


rdN 1 


dt + Mi dt 


■) 


dt 

-h t ^h 


-) t H 


Zi XrfAdZi ZjX d B QM ZfXdBi'd 
C 0 , c ,Z 1 Doom D 0hc , 

Ci t ciZi Dio, cl Du d 


' SA 

0 

SB 0 

SBi ' 

H T A 

h t arh 

H t B 0 

H t Bi 

c 0 

C 0 RH 

Doo 

Doi 

Cl 

CiRH 

Dio 

Dn 

\ In 

o 1 




+ 


where 


0 B 


'p 2 


0 


0 D 


P 12 J 


C n 


B a 

Dn 



0 R 

0 

0 1 


O 

o, 

I 

o 


r a 

■c^aw 

Baw 


SARH 

O' 

^aw 

Baw_ 


0 

0 


+ 


’N SB 2 ' 


a n 

-c^-aw LJ aw 


\m t h ol 

.0 [0 4 


Caw D aw - 


1 

o 


Multiply diag {Zj , I, I, /} from the left side and its 
conjugate transpose from the right side of eq. (A3) 
and we get 


Jll 

J T 

-J-21 

Lf ' 

J 21 

J 22 

Ll 

. L l 

L 2 

-n 


< 0 


(AS) 


with the shorthand notation 

Jn = SA + A S + y pi — — 
^ dpt 


i=l 


dS 


.dX 1 


■ hi = H 1 A + Ai w + //' (R— + M-) 


dt 


dR N 


Joo = H t (AR + RA t - yp~)H 

ti d P* 


H - Bp^C a w 

T o _1_ dT 


AT nT 


u - 


Lo = 


B£ S + n; iw + WKCo 
BjS 
Cl 


(S 0 T + WKCoR)E + D T aw B^ + WKCa 


Bi H 


C\ R H -f Dpi 2 C a w 
By Schur complement, it is equivalent to 
Jii + £,'11 '/. 


Jn +L{B- l Li 

21 T u 2 


joi +^ll 


■J‘22 + l 2 n 1 l 2 


< o 


(A6) 


Clearly, the lower (3x3) matrix of inequality (A5) 
must be negative definite; this determines the feasible 
D a w Let the (2, 1) element be equal to zeros and 
one can solve for A aw . This also leads to decoupled 
LMIs from the inequality (A6). Then B„ w ,C aw terms 
can be solved from the (1, 1) and (2, 2) elements of the 
decoupled inequality (A6). Note that both inequalities 
have regular solutions. 27 

The (1,1) element of the above matrix inequality 
corresponds to LMI (9) after elimination of the 
variables B aw and D aw . It can also be shown that the 
(2,2) element is equivalent to LMI (8) by eliminating 
Caw: Daw- Q.E.D. 
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